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ABSTRACT 

We present a new method for deriving cosmological constraints based on the evolution of the baryon mass 
function of galaxy clusters, and implement it using 17 distant clusters from our 160 deg 2 ROSAT survey. The 
method uses the cluster baryon mass as a proxy for the total mass, thereby avoiding the large uncertainties of 
the M tot -T or M tot -Lx relations used in all previous studies. Instead, we rely on a well-founded assumption 
that the Mt,/M tot ratio is a universal quantity, which should result in a much smaller systematic uncertainty. 
Taking advantage of direct and accurate Chandra measurements of the gas masses for distant clusters, we find 
strong evolution of the baryon mass function between z > 0.4 and the present. The observed evolution defines a 
narrowband in the £l m -A plane, f2,„+0.23A = 0.41 ±0.10 at 68% confidence, which intersects with constraints 
from the Cosmic Microwave Background and supernovae la near Q, m = 0.3 and A = 0.7. 
Subject headings: cosmological parameters — galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

The growth of large scale structure in the near cosmologi- 
cal past is very sensitive to the value of the density parameter, 
fi m , and, weakly, to the cosmological constant, A (see, among 
others, Sunyaev 1971, Peebles 1980, Carroll, Press & Turner 
1992). One of the strongest manifestations of the growth 
of density perturbations is the formation of virialized objects 
such as clusters of galaxies (Press & Schechter 1974). There- 
fore, the evolution of the cluster mass function provides a sen- 
sitive cosmological test (e.g. Evrard 1989 and later works). 

Present-day theory provides machinery for accurate pre- 
dictions for the cluster mass function at any redshift, in any 
cosmology (see a brief review in §0}. Most of the difficul- 
ties with using cluster evolution as a cosmological probe are 
on the observational side. The foremost problem is the ne- 
cessity to relate cluster masses with other observables, such 
as the X-ray temperature. Even for the few brightest, best- 
studied clusters, the total mass within the virial radius — the 
quantity generally used in theory — is measured presently 
with large uncertainties (see, e.g., Markevitch & Vikhlinin 
1997 and Fischer & Tyson 1997 for representative uncertainty 
estimates for the X-ray hydrostatic and weak lensing meth- 
ods, respectively). To measure virial masses in a large, com- 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cam- 
bridge, MA 02138; avikhlinin@cfa.harvard.edu 
- Space Research Institute, Moscow, Russia 

3 European Southern Observatory, Karl-Schwarzschild-StraBe 2, D-85748 
Garching, Germany 

4 Dpto. de Astronomia y Astrofisica, Pontificia Universidad Catolica, 
Casilla 306, 22 Santiago, Chile 

5 Department of Physics and Astronomy, Ohio University, Athens, OH 
45701 

6 Istituto di Radioastronomia CNR, Via P. Gobetti, 101, 1-40129, Bologna, 
Italy 

Danish Space Research Institute, Juliane Maries Vej 30, 2100 Copen- 
hagen O, Denmark 

8 University of Hawaii, Institute for Astronomy, 2680 Woodlawn Drive, 
Honolulu, HI 96822 



plete sample of clusters is presently impractical. Therefore, 
one usually resorts to other, more easily measurable quanti- 
ties, which exhibit scaling relations with the total mass. The 
most widely used observable is the temperature of the intra- 
cluster medium. It is expected that T is related to the virial 
mass as Ma T 3 / 2 , and the available observations generally 
support this law (Nevalainen, Markevitch & Forman 2000; 
Horner, Mushotzky & Scharf 1999; Finoguenov, Reiprich & 
Bohringer 2001). To normalize this relation, virial masses 
still have to be measured in a representative sample of clus- 
ters. Any systematic uncertainties in the cluster mass mea- 
surements translate into the same uncertainty in the normal- 
ization of the M-T relation. 

The X-ray temperature functions (XTF) for both nearby and 
distant clusters have been used as a proxy for the mass func- 
tion by a number of authors, using the normalization of the 
M-T relation from different sources. An analysis of a sta- 
tistically complete local XTF was first performed by Henry 
& Arnaud (1991). The analysis of the local XTF yields the 
normalization (erg), in combination with f2 m , and slope of the 
power spectrum of density fluctuations on cluster scales. This 
type of study has since been applied to ever improving obser- 
vational data (see Evrard et al. 2002 for a recent summary). 

A first cosmological measurement using the evolution of 
the XTF at redshifts greater than zero is due to Henry (1997), 
who derived fi m w 0.5 ±0.15. A similar analysis was per- 
formed by Eke et al. (1998), Donahue & Voit (1999), Henry 
(2000), Blanchard et al. (2000), with updates on both obser- 
vational and theoretical sides as well as with different nor- 
malizations of the M-T relation. We also note the studies 
by Reichart et al. (1999) and Borgani et al. (2001) who mod- 
eled the X-ray luminosity functions for distant clusters using 
an empirical L-T relation and the usual theory for the XTF 
evolution. 

Theoretical prediction for the number density of clusters 
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with a given temperature is exponentially sensitive to the nor- 
malization of the M-T relation. A 30% change in the nor- 
malization of the M-T relation results in a 20% change in 
the determination of the erg parameter and shifts an estimate 
of fi m by a factor of 1.3-1.4 (e.g., Rosati, Borgani & Norman 
2002). Given that ±50% is a fair estimate of the present sys- 
tematic uncertainties in the M-T normalization, the impact 
of this relation on cosmological parameter determination can 
be very significant. 

In modeling the high-redshift temperature function, the the- 
oretically expected evolution of the M-T relation is usually 
assumed: M oc r 3 / 2 / (1 + z) 3//2 . It would be very important to 
verify this evolution at high z. Unfortunately, this seems to 
push the limits of current observational techniques. The evo- 
lution of other scaling relations can be also quite important. 
For example, the L—T relation is required for computing the 
volume of X-ray flux limited surveys, and for estimating clus- 
ter masses, when the X-ray luminosity function is used as a 
proxy for the mass function. Most of the previous studies 
have assumed that the L— T relation does not evolve, however 
recent Chandra measurements indicate otherwise (Vikhlinin 
et al. 2002). 

Finally, most of the available cosmological studies are 
based on the distant cluster sample from the Einstein Ex- 
tended Medium Sensitivity Survey (Gioia et al. 1990, Henry 
et al. 1992). The completeness of this survey at high redshifts 
is often questioned (Nichol et al. 1997; Ebeling et al. 2000; 
Lewis et al. 2002). Although we believe that the EMSS is es- 
sentially complete (as confirmed by good agreement between 
the XLFs from the EMSS and several ROSAT surveys — 
Gioia et al. 2001, Vikhlinin et al. 2000, Mullis et al., in prepa- 
ration), it is arguably important to try another completely in- 
dependent sample for the cosmological studies based on clus- 
ter evolution. 

The purpose of this paper is to apply the evolutionary test 
to a new sample of distant, z > 0.4, clusters derived from the 
160 deg 2 ROSAT survey (Vikhlinin et al. 1998), relying on 
the evolution of the cluster scaling relations as measured by 
Chandra, and using a modeling technique which does not rely 
upon the total mass measurements to normalize the M-T or 
M-L relations and therefore bypasses many of the uncertain- 
ties mentioned above. 

The proposed method relies on the assumption that the 
baryon fraction within the virial radius in clusters should be 
close to the average value in the Universe, fay = Sit /^m- This 
is expected on general theoretical grounds (White et al. 1993) 
because gravity is the dominant force on cluster scales. The 
universality of the baryon fraction in massive clusters is sup- 
ported by cosmological numerical simulations (e.g. Bialek, 
Evrard & Mohr 2001) and by most available observations 
(most recently, Allen et al. 2002). In principle, the equality 
fi, = fib/ ft m can be used to determine cosmological parame- 
ters either from the absolute measurements of fb in the nearby 
clusters (White et al. 1993) or from the apparent redshift de- 
pendence of fb measurements in the high-redshift clusters 
(Sasaki 1996, Pen 1997). Both methods, however, require 



an accurate observational determination of the total mass in 
individual clusters, which we would like to avoid. 

Unlike the total mass, the baryon mass in clusters is rel- 
atively easily measured to the virial radius (Vikhlinin et al. 
1999). To first order, the baryon and total mass are triv- 
ially related, Mb = MQ,b/Sl m . If this is the case, the relation 
between the cumulative total mass function, F(M), and the 
baryon mass function, Fb(Mb), is also trivial, 

Fb{M b ) = F(n m /n b Mb). (i) 

The average baryon density in the Universe is accurately 
given by observations of the light element abundances and Big 
Bang Nucleosynthesis theory (Buries, Nollet & Turner 2001). 
Therefore, we have everything needed to convert the theoreti- 
cal model for the total mass function, F(M), to the prediction 
for a directly measurable baryon mass function — to com- 
pute F(M) one must choose a value of 57,„, which fixes the 
scaling between the total and baryon mass. The conversion 
of the total to the baryon mass function — at least to a first 
approximation — can be expressed through the model param- 
eters, and therefore Fb(Mb) can be modeled with the usual set 
of parameters, fl m , a%, n, h, and at a high redshift — A. We 
will show below that the evolution of the baryon mass func- 
tion provides robust constraints on a combination of il m and 
A. Modeling of the baryon mass function at z = provides 
a measurement of the power spectrum normalization, <j%, and 
of the shape parameter T = fl m h, with degeneracies that differ 
from those given by the temperature function (Voevodkin & 
Vikhlinin, in preparation). 

Theoretical models of the cluster mass function usually 
consider masses measured within radii defined by certain val- 
ues of the density contrast, S = 3M(r)/(47rr 3 )//?, where p is 
either the critical density or mean density of the Universe. The 
assumption that the baryon fraction in clusters is universal at 
sufficiently large radii allows one to determine the overdensity 
radii r$ without measuring the total mass. Indeed, the baryon 
and matter overdensities are then equal, 8b = S, and the baryon 
overdensity is defined by a baryon mass Mb{r) and the mean 
baryon density in the Universe, known from the BBN. One 
of the most accurate theoretical models for the mass function 
(Jenkins et al. 2001) uses masses corresponding to the over- 
density 5 = 324 relative to the mean density of the Universe 
at the redshift of observation. The matching baryon masses, 
Mfc 324 are easily measured using ROSAT data for nearby clus- 
ters (Voevodkin, Vikhlinin & Pavlinsky 2002a, VVPa here- 
after), and Chandra data for distant clusters (Vikhlinin et al. 
2002). 

The paper is organized as follows. In §|2] we discuss the 
distant cluster sample and the gas mass measurements in the 
distant and nearby clusters. The baryon mass function for dis- 
tant clusters is derived in §[3] The theory for modeling these 
data is reviewed in §|4] We describe our fitting procedure in 
§0 and present the derived constraints on £1„, and A in §[6] 

All cluster parameters are quoted for h = 0.65, f2 m = 0.3, 
and A = 0.7. The X-ray fluxes and luminosities are in the 0.5- 
2 keV energy band. Measurement uncertainties are la. 
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2. DATA 

2.1. Distant Cluster Sample 

The sample of distant clusters used in this work is derived 
from our 160 deg 2 survey (Vikhlinin et al. 1998). This sur- 
vey is based on the X-ray detection of serendipitous clusters 
in the central parts of a large number of the ROSAT PSPC 
observations. The spatial extent of the X-ray sources is used 
as a primary selection criterion, but it is verified by indepen- 
dent surveys (Perlman et al. 2002) that above the advertised 
flux limits, we do not miss any clusters because of the lim- 
ited angular resolution of ROSAT. Therefore, for all practical 
purposes the 160 deg 2 cluster sample can be considered as 
flux-limited. 

The optical identification of the entire X-ray sample is now 
complete and redshifts of all clusters are spectroscopically 
measured (Mullis et al., in preparation). Optical identifica- 
tions show that our X-ray selection is very efficient. In the 
entire sample, only 20 of the 223 detected X-ray cluster can- 
didates lack a cluster identification. Above the median flux of 
the survey, / > 1.4 x 10~ 13 erg s" 1 cm" 2 , 111 of 114 sources 
(or 98%) were confirmed as clusters. Therefore, the 160 deg 2 
survey is in effect purely X-ray selected. 

Our X-ray detection procedure is fully automatic, and there- 
fore all essential statistical characteristics of the survey can be 
derived from Monte-Carlo simulations. We have performed 
very extensive simulations to calibrate the process of detec- 
tion of clusters with different fluxes and radii. This gave an 
accurate determination of the survey sky coverage as a func- 
tion of flux (Vikhlinin et al. 1998), and therefore of the survey 
volume for any combination of cosmological parameters. 

2.2. Chandra Observations 

As of Fall 2002, Chandra has observed 6 clusters from 
the 160 deg 2 sample. These objects almost complete a flux- 
limited subsample, / > 2 x 10~ 13 erg s" 1 cm" 2 , of objects at 
z > 0.4 (7 total), and therefore they are suitable for statistical 
studies. Chandra observations of these clusters are discussed 
in Vikhlinin et al. (2002). Briefly, the quality of the X-ray data 
allows a measurement of the cluster temperature to a 10% ac- 
curacy, and of the gas mass at the virial radii with a 20% un- 
certainty. For the purposes of this paper, we need only the 
measurements of the gas mass. 

The Chandra observations of the distant 160 deg 2 clusters 
as well as those from the EMSS and RDCS (Rosati et al. 
1998) samples show that scaling relations between the X-ray 
luminosity, temperature, and the gas mass significantly evolve 
with redshift (Vikhlinin et al. 2002). In particular, the L-M g 
relation evolves so that the gas mass for a fixed luminosity is 
proportional to (1 + z) _1 ■ 8±0 - 25 but the slope and scatter in the 
high-redshift relation are the same as in the local relation ( Vo- 
evodkin, Vikhlinin & Pavlinsky 2002b, VVPb hereafter). This 
allows an estimate of the gas mass, with a ±35% uncertainty, 
from the easily measured X-ray flux. 

The L-M g relation allows us to extend the estimate of the 
gas mass to all high-redshift clusters that do not have Chan- 
dra observations. However, we will use only clusters with 



f > 1 .4 x 10" 13 erg s" 1 cm" 2 . At z = 0.4, this flux limit corre- 
sponds roughly to T = 3 keV. It is possible that for much less 
massive clusters, our basic assumption of the universality of 
the baryon fraction may not hold due to, e. g., preheating of 
the intracluster medium (e.g., Bialek et al. 2001). The median 
redshift for the selected distant clusters is z = 0.55. 

The measurements and estimates of the baryon mass in dis- 
tant clusters are listed in Tabled The measured masses were 
derived using direct deprojection of the X-ray surface bright- 
ness profiles. The mass uncertainties were obtained by error 
propagation in the process of deprojection and so correctly re- 
flect the statistical noise in the surface brightness profiles near 
r324. The only object with Chandra data which was not used 
by Vikhlinin et al. (2002), is CL 0030+2618. This cluster was 
observed very soon after the Chandra launch while the detec- 
tor parameters did not yet reach their nominal values. These 
data are unusable for spectral (temperature) measurements, 
but the imaging analysis which yields the gas mass can be 
performed without any complications. 

The baryon masses were estimated from the M/, -L correla- 
tion for approximately 60% of clusters in our distant sample. 
Both quantities in the Mh - L correlation are easily and di- 
rectly measured, and so there is no systematic uncertainty in 
the derivation of normalization, or slope, or evolution of this 
relation. 

Finally, we remark on the excellent agreement of the X-ray 
fluxes for the 160 deg 2 survey clusters derived from ROSAT 
and Chandra data. For the 6 clusters in common, the ROSAT 
fluxes deviate by no more than 13%, always within the sta- 
tistical uncertainties. The average ratio of the ROSAT and 
Chandra fluxes is 0.99 ± 0.06. 

2.3. Summary of the Low-Redshift Results 

The baryon mass measurements in a large sample of low- 
redshift clusters, as well as the determination of the baryon 
mass function at z = 0.05 was reported in VVPa,b. These re- 
sults provide a low-redshift reference point for our measure- 
ment of cluster evolution, so a brief summary is in order. 

VVP have used a flux-limited subsample of the low-redshift 
clusters detected in the ROSAT All-Sky Survey. The initial 
object selection was performed using the HIFLUGCS sam- 
ple (Reiprich & Bohringer 2001). The X-ray fluxes were re- 
measured and 52 clusters were selected in the redshift interval 
0.01 < z < 0.25 and with f x > 1.4 x 10" 11 erg s" 1 cm" 2 . The 
median redshift of this sample is z = 0.05. Most of these clus- 
ters have ROSAT PSPC pointed observations which were used 
to measure the gas mass (the procedure is described in VVPa 
and Vikhlinin, Forman & Jones 1999). 

Using a subsample of clusters with published optical mea- 
surements, a tight correlation between the gas mass and op- 
tical luminosity of the cluster has been established. This re- 
lation allows an estimate of the baryon (i. e. gas+galaxies but 
excluding intergalactic stars) mass of all clusters from the X- 
ray data alone: 

-0.5 

1.100 + 0.045 



M h .324 = M„ 32 4 X 



g,324 



10 15 M, 
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Table 1 

x-ray data and baryon masses for the distant 160 deg 2 clusters 



Zo.5-2 (ROSAT) Lo.5_2 M g i2 A *4 i32 4 /0.5-2 (Chandra) 

Cluster z (ergs' 1 cm' 2 ) (ergs' 1 ) (10' 3 M e ) (1O 13 M ) (ergs' 1 ) 

CL 1416+4446 0.400 4.04 x 10' 13 2.5 x 10 44 8.0 ±2.3 9.8±2.8 4.51 x 10' 13 

CL1701+6414 0.453 3.86 x 10' 13 2.9 x 10 44 10.4 ±2.8 12.5 ±3.4 4.14xl0' 13 

CL1524+0957 0.516 3.04 x 10' 13 2.7 x 10 44 9.3 ±2.3 11.2±2.8 2.81 x 10' 13 

CL 1641+4001 0.464 2.94 x 10' 13 2.3 x 10 44 ■■■ 8.3 ±2.9 

CL0030+2618 0.500 2.43 x 10' 13 2.3 x 10 44 6.6 ±2.7 8.1 ±3.3 2.51 x 10' 13 

CL1120+2326 0.562 2.13 x 10' 13 2.2 x 10 44 6.6± 1.5 8.2± 1.9 1.88xlO' 13 

CL1221+4918 0.700 2.06 x 10' 13 4.1 x 10 44 11.3 ±2.1 13.6±2.6 2.25 x 10' 13 

CL 0853+5759 0.475 1.98 x 10' 13 1.6 x 10 44 ~ 6.3 ±2.2 ~ 

CL0522-3625 0.472 1.84xl0' 13 1.5 x 10 44 ■■■ 5.9±2.1 

CL1500+2244 0.450 1.78xl0' 13 1.3 x 10 44 ■■■ 5.5±1.9 

CL0521-2530 0.581 1.76 xl0' 13 2.3 x 10 44 ■■■ 7.2±2.5 

CL0926+1242 0.489 1.67xl0' 13 1.5 x 10 44 ■■■ 5.7±2.0 

CL0956+4107 0.587 1.56 xlO' 13 2.0 x 10 44 ■■■ 6.6±2.3 

CL1216+2633 0.428 1.54xl0' 13 l.OxlO 44 ■■■ 4.6±1.6 

CL1354-0221 0.546 1.45 X 10' 13 1.6xl0 44 ■■■ 5.8±2.0 

CL1117+1744 0.548 1.44 xl0' 13 1.6xl0 44 ■■■ 5.8±2.0 

CL1213+0253 0.409 1.43 x 10' 13 0.9 x 10 44 ■■■ 4.2 ±1.5 



All distance-dependent quantities are derived assuming Q m = 0.3, A = 0.7, and h = 0.65. For different values of the Hubble constant, masses scale as h . 
The luminosities are determined using the ROSAT fluxes, not the Chandra measurements. The baryon mass (i. e. gas+stars) is estimated from the gas mass 
using J2}. For CL 1641+4001 and all clusters with / < 2 X 10~ 13 erg s' 1 cm' 2 , the values of 324 are estimated from the L—Mf, relation. 



where Mb.yiA is the baryon mass corresponding to the over- 
density S = 324, and M g ^iA is the gas mass corresponding to 

5 = 324 without accounting for the stellar mass. The stellar 
contribution to the baryon mass is small, but non-negligible, 
of order 10-15% for massive clusters (this is consistent with 
the estimates in Fukugita, Hogan & Peebles 1998). This is es- 
timated with approximately 25% uncertainty (i. e. ±3% of the 
total baryon mass). Assuming that the ratio of stellar and gas 
mass does not evolve at high z, i. e. that the galaxies neither 
confine the intracluster gas nor lose their mass significantly 
due to stellar winds, equation © can be used to estimate the 
baryon mass in distant clusters (our results are not sensitive to 
this assumption). 

VVPb have derived the low-redshift baryon mass function 
(reproduced below). The measurement uncertainties in this 
mass function were determined using a technique which ac- 
counts for the Poisson noise as well as the individual Mb mea- 
surement uncertainties (see also Appendix). 

To summarize, the baryon mass function for the low- 
redshift clusters is measured reliably. By itself, it provides 
constraints on the shape and normalization of the power spec- 
trum of the density perturbations, i.e. on the parameters a% 
and r = VL m h. These results will be reported in Voevodkin 

6 Vikhlinin (in preparation). In the present Paper, we re- 
strict ourselves to the cosmological constraints derived from 
the evolution of the mass function between z > 0.4 and z — 0. 

3. BARYON MASS FUNCTION AT Z > 0.4 

In this section, we use the Chandra observations to derive 
the baryon mass function for distant clusters. First, we de- 
scribe the computations of the survey volume as a function of 



3.1. Survey Volume 

We consider the mass function for 0.4 < z < 0.8. Since our 
cluster sample is selected by X-ray flux, the computation of 
the surveyed volume as a function of mass is only possible 
if there is a correlation between the baryon mass and X-ray 
luminosity. Such a correlation is observed (VVPb, Vikhlinin 
et al. 2002), 

M,,cx {1+z) Aml L x * 3 (3) 

with a 42% scatter in luminosity and 35% scatter in mass; Lx 
is the total luminosity in the 0.5-2 keV band. 

In the case of zero scatter in the Mb~L relation, the volume 
is computed as 

(1+z)- l20Aml S(z) 



V = 



dz 4ndi{z) 



(4) 



where A(f) is the survey area as a function of flux, dV /dz is 
the cosmological dependence of the comoving volume, equa- 
tion {3} is rewritten as L = CM L20 (1 + z)~ l20Aml , and S(z) is 
a correction for redshifting of the bandpass (equivalent to the 
K-correction in the optical astronomy). Note that for the non- 
evolving L — T relation, S(z) is independent of cosmology. 
The corrections due to the observed rate of evolution in this 
relation (Vikhlinin et al. 2002) are small, ~ 5% in volume, 
and they were ignored. 

Introducing the scatter in the M^-L relation with the log- 
normal distribution (see VVPb), equation @ becomes 

(^£'£ exp( " is3r/2< ' 2) '' igrx _ 

/ CM L20 {l+z)- L2{)A ^S(z)\ dV D) 
x A [y — 1 



V = 



4ndi(z) J dz dZl 

where a = 0.173 is the observed rms log-scatter in the local 
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FIG. 1 . — Left: Comoving volume of the 160 deg 2 survey for 0.4 < z < 0.8 as a function of the cluster baryon mass. Right: Volume per unit redshift interval 
for three values of Mt, — 0.25, 0.5, and 1 X 1O I4 M0 (from bottom to top). On each panel, the solid lines are for the fl m = 0.3, A = 0.7 cosmology, and dashed 
lines are for Q„, = 1 , A = 0. The volumes are computed for a limiting X-ray flux of / m ; n = 




mean Mb~L relation (VVPb). Numerical integration of this 
equation gives the survey volume as a function of baryon mass 
for any set of the cosmological parameters fi m , A, and h. The 
inner integral in l|5} represents the effective value of the survey 
volume per unit redshift interval for clusters of given mass, 
(dV/dzU- 

There is no simple scaling of V(M) with either f2 m or 
A. In addition to dV /dz, the index A ML is also cosmology- 
dependent. When fitting the mass function to different cos- 
mological models, we recomputed V (M) for every new com- 
bination of the model parameters. 

Figure[2shows the comoving volume in the redshift interval 
0.4 < z < 0.8 as a function of the baryon mass. Note that in 
the limit of large masses, the survey volume approaches that 
of the local all-sky surveys. 

The right panel of Fig^shows a volume per unit redshift 
interval for three values of Mb- The clusters with small masses 
generally have low luminosities and therefore can be detected 
only near the lower boundary of the redshift interval consid- 
ered. Massive clusters can be detected at any z, and there- 
fore most of the volume for such clusters is accumulated near 
the higher boundary. This dependence of (dV/dz) e f£ on mass 
should be taken into account in modeling the mass function 
because the evolutionary effects between z = 0.4 and 0.8 are 
quite strong. 



3.2. Observed Baryon Mass Function 

Since the number of clusters in our sample is rather small, 
we will work with the integral representation of the mass 



function, which is determined as 
F{M) =N(> M) = 



E 

M, (,) >M 



V{Mf 



The Poisson errorbars for F(M) are computed as 
(A^) 2 = E ' 



(6) 



(7) 



Ml ' >M 

In addition to the Poisson noise, the F(M) errorbars should in- 
clude an additional component reflecting the individual mass 
measurement uncertainties. This was done using the Max- 
imum Likelihood analysis described in the Appendix. We 
estimate that other obvious sources of uncertainty in F (M) 
are small. For example, the inaccuracies of the survey vol- 
ume caused by uncertainties in the 160 deg 2 sky coverage 
or variations of the evolution or scatter in the Mb - L rela- 
tion within the observationally allowed intervals are all below 
10%, which is much smaller than the Poisson noise in the 
mass function, so we ignore them. 

The baryon mass function computed in the Sl„, = 0.3, A = 
0.7 cosmology is shown in Fig. |2] A comparison with the 
local F (M) from VVPb shows a strong evolution — at z « 0.5 
there are approximately 10 times fewer clusters of the same 
mass per unit comoving volume than at z = 0. This evolution 
is quite insensitive to the assumed cosmology. The right panel 
in Fig.|2]shows the mass functions computed assuming several 
combinations of £l m and A. Mostly, the changes are restricted 
to a translation along the band defined by F(M) itself. 

Part of the evolution of the mass function in Fig.|2]may be 
unrelated to the physical growth of clusters and instead re- 
lated to the commonly used definition of mass. Indeed, let us 
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Fig. 2. — Left: The baryon mass function for the 160 deg 2 survey sample in the redshift interval 0.4 < z. < 0.8. The local mass function from Voevodkin et al. 
(2002b) is shown by triangles. The grey band shows a 68% uncertainty interval for the mass functions, including both the Poisson noise and mass measurement 
uncertainties. The mass function at z. > 0.4 is shown by circles. The solid circles correspond to clusters with Chandra baryon mass measurements, and the 
empty circles to the mass estimates from the Mj —L relation. The high-redshift mass function is computed for £l m = 0.3, A = 0.7; the local mass function is 
cosmology-independent. Right: Baryon mass functions derived assuming f2 m = 0.3, A = 0.7 (empty circles); O m = 0.3, A = (solid circles); Q,„ = 1, A = 
(squares). The solid line shows a fit to the local mass function. On both panels, the data are for h = 0.5. The measured M^, scales as ex. /j~ 2 ' 25 and volume 
scales as usual, V oc /j~ 3 



assume that all clusters stay exactly the same at all redshifts. 
The operational definition of mass used in our study is related 
to the density contrast above the mean density of the Uni- 
verse at the cluster redshift. Since the mean density changes 
as (1 +z) 3 , the radius of the fixed contrast decreases with red- 
shift, and therefore the mass also decreases. For the average 
density profile at large radii, we find that this effect results in 
the mass scaling M ~ (1 +z)~ 1,5 . This mass offset explains 
about 40% of the observed effect, and the rest is due to the 
real growth, by a factor of ~ 2 in mass, of clusters between 
z > 0.4 and z = 0. 

As was explained above, the total mass functions are re- 
lated to the baryon mass functions through a constant log shift 
along the mass axis, which does not alter the vertical offset 
between the local and high-redshift functions. The vertical 
offset is related to the growth of density fluctuations between 
z > 0.4 and z = 0. Below, we use it to derive constraints on the 
parameters il m and A. 



4. THEORY 

Our computations of the theoretical mass functions are 
based on the Jenkins et al. (2001) universal fit to the mass 
functions in large-volume cosmological simulations. Jenkins 
et al. find that the mass functions expressed in terms of the 
rms density fluctuations on mass scale M are of a universal 
form: 

M dn(M) 



/W=-77v 



(p) din a 



= Aexp(--|ln<7 +b\ c ) 



(8) 



where n(M) is the number of clusters with mass greater than 
M. The universality means that the coefficients A, b, c do not 
depend on the cosmological parameters, nor on redshift. All 
the cosmological dependencies are in the function a(M,z). 
The values of A, b, and c depend, however, on the definition 
of the cluster mass. For the definition in use here — mass 
corresponding to the spherical overdensity 5 = 324 relative to 
the mean density of the Universe at the redshift of observation 
— Jenkins et al. provide the values 

A = 0.316, £ = 0.67, c = 3.82. (9) 

The function a(M,z) is the product of the present-day rms 
density fluctuations, a(M), and the growth factor of linear 
density perturbations, D(z). The growth factor as a function of 
z is uniquely determined by the values of the cosmological pa- 
rameters n m and A. We used the software provided by Hamil- 
ton (2001) to compute D(z) in any cosmology. The function 
<t(M) is trivially computed from the present-day linear power 
spectrum, P(k). We assume that P(k) is the product of the 
inflationary power law spectrum, k", and the transfer func- 
tion for the given mixture of CDM and baryons. The com- 
putation of the transfer function for any set of cosmological 
parameters was performed using the analytic approximations 
developed by Eisenstein & Hu (1998) to the exact numerical 
calculations with CMBFAST (Seljak & Zaldarriaga 1996). 

Finally, the computation of the theoretical mass function ac- 
cording to (|8} was performed using the software kindly pro- 
vided by A. Jenkins. The program was slightly modified to 
include the power spectrum model by Eisenstein & Hu and 
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the D(z) computations due to Hamilton. With these modifica- 
tions, we have the machinery for precise computation of the 
total cluster mass function for any combination of the cosmo- 
logical parameters at any redshift. 

As was discussed above, if the baryon fraction in clusters 
equals its universal value, the model baryon and total mass 
functions are related as Fi,(Mb) = F(Mbfl m /fli,)- Let us con- 
sider how this is modified if the baryons in clusters are under- 
abundant by a factor T: Mb/M = Tilb/fl m . In this case, 
the total mass that corresponds to the measured baryon mass 
Mb,324, equals M324 = Y~ L5 J7 m /f4A^,324- The factor T" 1 here 
reflects the underabundance of baryons. The additional factor 
T"° 5 follows from the fact that the total and baryonic over- 
densities are now different. We would like to have the masses 
within the radius of total overdensity 5 = 324. This corre- 
sponds to the baryonic overdensity 5b = T 324, but we mea- 
sure Mb within 5/, = 324. We find empirically that the den- 
sity profile of a typical cluster corresponds to the mass profile 
Mb ,5 oc (5" ' 5 , and hence the measured baryon masses have to 
be scaled approximately by the factor T" 5 . Therefore, we 
have the following relation for the theoretical model of the 
experimental baryon mass function: 

F h (Mb)=F(T- L5 n m /n h Mb), (10) 
where F(M) is computed from equation 0. 

4.1. Dependence of the theoretical mass function on the 
model parameters 

Full details of our modeling procedure will be given below 
(§0. Here we present only a brief outline necessary to under- 
stand the influence of the model parameters on the proposed 
cosmological test. Our measurements are the baryon mass 
functions at z ~ and z = 0.4 — 0.8. All the constraints will 
be obtained by the requirement that the model mass functions 
are consistent with both the low- and high-redshift measure- 
ments. In this Paper, we consider only constraints from the 
cluster evolution and minimize the usage of any other infor- 
mation. In particular, we will not fit the shape of the local 
baryon mass function, but only use its normalization in the 
5 x 10 13 M Q <Mb < 10 14 M Q mass range (approximately the 
median mass of the low-redshift sample). The above mass 
interval corresponds to ICM temperatures of 3.7-5.5 keV. 

The theoretical model for the baryon mass function depends 
on the following parameters: il m ; A; average baryon density, 
Flbh 2 ; the slope of the primordial power spectrum, n; present- 
day normalization of the power spectrum, a%; Hubble con- 
stant, h; and also on any deviations of the baryon fraction in 
clusters from the universal value. Let us discuss how each of 
these parameters enters the model computations and compar- 
ison of the model with observations. 

Parameter n: The slope of the primordial power spectrum 
determines the general slope of the mass functions. When n 
varies, the high-redshift and local mass functions change in 
unison. Since in the nearby clusters, we consider only the 
overall normalization of the mass function, and the measure- 
ment of the high-z mass function has large statistical uncer- 
tainties, our evolutionary test is insensitive to variations of n 



in the interval 0.5 < n < 1 .5 (as we explicitly verified). This is 
much wider than the interval allowed by the current CMB ob- 
servations, n = 1.05 ±0.1 (Sievers et al. 2002, Rubino-Martin 
et al. 2002). We use n = 1 for our baseline model. 

Parameter h: The Hubble constant enters through its ef- 
fect on the shape of the power spectrum (mainly via the prod- 
uct il m h, Bond & Efstathiou 1984). It also affects the scal- 
ing between the total masses in the model and the measured 
baryon masses of clusters. The effect of h on the power spec- 
trum is unimportant for our test for the same reason, as that for 
n. The effect of h on the scaling between Mb and M is a shift 
of the model mass function by a factor of h 015 with respect 
to the measurements (for details, see Voevodkin & Vikhlinin 
in preparation). In the mass range of interest, the model mass 
functions at high z and at z = are almost parallel and there- 
fore a shift along the mass axis can be compensated by a slight 
change in a% . We have checked, that the constraints on fl m and 
A are unchanged for any h in the interval 0.5 < h < 0.8. For 
the baseline model, we use h = 0.65. 

Parameter flhh 2 : The average baryon density in the Uni- 
verse is given by BBN theory, n b h 2 = 0.020 ±0.001 (Buries et 
al. 2001). A larger variation of this parameter would cause a 
slight change in the shape of the present-day power spectrum, 
and also would affect the scaling between the total and baryon 
mass. As was already noted, these changes are unimportant 
for our evolutionary test. We have verified that our results 
are unchanged for any value in the interval 0.015 < ilbh 2 < 
0.025. In the baseline model, we use VLi,h 2 = 0.02. 

Parameter a%: The normalization of the mass functions 
is exponentially sensitive o%. The requirement that the the- 
oretical baryon mass function fits the measurements at z = 
effectively fixes cr% for any combination of other model pa- 
rameters. The measurement of erg by this technique will be 
presented elsewhere. For the purposes of this study, all our 
constraints are marginalized over the acceptable (i.e., those 
consistent with the local mass function) values of erg. 

Parameters VL m and A: The cosmological density param- 
eter determines (mainly via the product il m h) the slope of the 
power spectrum on cluster scales, and hence the shape of the 
mass function. By design, our test is rather insensitive to the 
shape of the mass function (see discussion of the parameter n 
above). 

In addition, a combination of parameters il m and A deter- 
mines the growth factor of linear density perturbations, and 
therefore the evolution of the mass functions between high 
redshift and z = 0. This effect is the basis for our evolutionary 
test. 

Deviations of the baryon fraction from universality. Any 
deviations of the cluster baryon fraction from the universal 
value, flb/£l m , can potentially be a serious problem for cos- 
mological tests using the baryon mass function. There are 
three possibilities: 1) T/l, but is the same for all clusters; 2) 
T = 1 on average, but there is a cluster-to-cluster scatter; 3) T 
is a function of cluster mass (e. g. T becomes small for low- 
mass clusters). Taking these possibilities into account is im- 
portant for tests based on the shape of the baryon mass func- 
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tion and for the measurement of a% (a detailed discussion will 
be given in Voevodkin & Vikhlinin). However, the evolution- 
ary test discussed here is rather insensitive to T, as long as 
there is no strong evolution of this parameter with redshift, as 
shown below. 

Let T 7^1, but be the same for all clusters. This could be due 
to different virialization processes for the dark matter and the 
baryon gas. Adiabatic numerical simulations of cluster for- 
mation often lead to T = 0.9-0.95, independent of the cluster 
mass (Mathiesen, Mohr & Evrard 1999). That paper also con- 
siders the systematic bias in the X-ray measurements of the 
baryon mass due to deviations of cluster shapes from spheri- 
cal symmetry, which for our purposes is indistiguishable from 
T/l. The Mathiesen et al. results suggest that deviations 
from spherical symmetry are equivalent to T = 1.02- 1.12, 
almost independent of mass. The effect of a constant T ^ 1 
on our evolutionary test is to shift all mass functions along 
the M axis by T" L5 (cf. eq. II 01. As was already discussed, 
this shift is unimportant. We verified that the obtained cos- 
mological constraints are virtually identical for T = 0.85, 1, 
and 1.15. 

If there is scatter in T, the model for the baryon mass func- 
tion should be convolved with some kernel. As long as the 
scatter is not too large (the observational upper limits on the 
variations of fb are approximately 15-20%, see, e.g., Mohr et 
al. 1999), its effect on the mass function is very small. 

Let us now consider the possibility of a systematic change 
of T with the cluster mass. For example, a strong preheating 
of the intergalactic medium can cause the baryons to be un- 
derabundant in low-mass clusters (Cavaliere, Menci & Tozzi 
1997). Bialek et al. (2001) discuss the cosmological simula- 
tions of clusters with the level of preheating adjusted to re- 
produce the local L x - T relation. They find a systematic de- 
crease of T in low-temperature clusters, which can be written 
in terms of the baryon mass approximately as: 

T=l-M*/M fe , M* = 0.9x 10 13 M Q . (11) 

Such a dependence of T on mass distorts the shape of the 
experimental baryon mass function quite significantly, but 
mainly for the low mass clusters which we do not use here. 
We used eq. il It in our baseline model, but have verified that 
with T = 1, the results are virtually identical. 

Finally, we also considered a possible evolution of eq. dl It 
with redshift. A reasonable assumption is that M* in dl It 
corresponds to the fixed ratio of the entropy of the intra- 
cluster and intergalactic media, which leads to the scaling 
M» oc (1 +z) 3//2 . We verified that the constraints on fl m and 
A are very weakly affected by such a scaling. 

The evolution of T is important for our test only if it is sig- 
nificant for massive clusters. The amplitude of the statistical 
uncertainties on fi m and A is equivalent to a ±13% change in 
T at z = 0.55. We consider such changes unlikely because the 
physically motivated value for massive clusters is T = 1. 

5. FITTING PROCEDURE 

The tightest parameter constraints can be obtained by the 
maximum-likelihood analysis of the distribution of clusters as 



a function of both mass and redshift. However, given the nov- 
elty of our cluster evolution test, we decided to model the two 
experimental mass functions, at z ~ and 0.4 < z < 0.8, using 
goodness of fit information and rejecting only those models 
which are clearly inconsistent with the data. This is a simpler 
and more convincing method then the maximum-likelihood 
fit. Our approach is described below. 

While the computation of the model at z = is straightfor- 
ward (eq.|8]|9]^|[^at z = 0), it is less so for our high-z sam- 
ple because one expects a significant evolution within the red- 
shift interval considered, 0.4 < z < 0.8. To take this into ac- 
count, we subdivided this interval into 10 narrow z bins, com- 
puted cumulative mass functions on the same grid of masses 
in each bin, and then weighted them with (dV /dz) e ff (see eq.|5] 
and Fig. Q for each mass. This gives the model mass func- 
tion which can be directly compared with the one observed at 
0.4<z<0.8. 

Note that the volume computations depend on the evolution 
of the Mb~L relation (cf. eq.[5}, which is slightly cosmology- 
dependent. For self-consistency, for each combination of f2,„ 
and A we recomputed masses and luminosities of all clusters 
in the Vikhlinin et al. (2002) sample and rederived the value 
of the index Ami (cf. eq. [3J- This index also enters the vol- 
ume computations for the experimental baryon mass function, 
which we also recomputed for each set of parameters fi m , A, 
and h. 

As was discussed above, our test is insensitive to the precise 
values of the cosmological parameters n, fib, and h. There- 
fore, we fixed them at n = 1, h = 0.65, and Vlb = 0.02 h 2 , 
and varied only ft,,,, A, and erg. The model mass func- 
tions were computed on a grid fi m = 0.05,0.075, 1.25, 
A = 0.0,0.05,..., 1.1, and a% = 0.4, 0.405,..., 1.5. The con- 
straint on fl m and A was marginalized over a% — a combina- 
tion of Q, m and A was considered acceptable if for some value 
of erg, the model mass functions were consistent with the local 
and high-redshift measurements simultaneously. 

The comparison of the theoretical mass functions with ob- 
servations was designed to rely mainly on the observed evo- 
lution and to limit the use of any other information. In par- 
ticular, we use only the normalization of the local mass func- 
tion in the mass range 5 x 10 13 M Q < M b < 1 x 10 14 M Q (for 
h = 0.65). In the distant sample, we will consider the F(M) 
uncertainty intervals in the identical mass range. 

A simple and conservative comparison of the theoretical 
and observed mass functions is performed as follows. If the 
model cumulative function is entirely within the 68% confi- 
dence interval of the observed one, it is considered acceptable 
at the 68% confidence level. The confidence levels of 90% 
and 95% are applied identically. The combined confidence 
level for the local and distant mass functions are combined 
according to Table [5] For example, a model is acceptable at 
the 90% confidence level if it is within the 68% error bars of 
one of the mass functions, and within the 90% error bars of 
another. It can be shown that the described method leads to 
more conservative (wider) confidence intervals compared to 
the likelihood ratio test. 
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Table 2 

Combined confidence level for the local and 
high-redshift mass functions 



z = 

z>0.4\^ 


68% 


90% 


95% 


68% 


68% 


90% 


95% 


90% 


90% 


95% 




95% 


95% 







6. RESULTS 

The resulting 68% and 90% confidence regions for param- 
eters Q, m and A are shown in Fig. [3] The observed rate of the 
evolution of the baryon mass function defines a narrow band 
in the (f2 m , A) plane, slightly inclined to the A axis: 

fi m + 0.23A = 0.41 ±0.10. (12) 

The uncertainty in Q m for any fixed A is ±0.10 at the 68% 
confidence level. For a flat cosmology, tt m + A = 1, we have 
fl m = 0.24 ±0.12. Using only clusters with direct Chandra 
measurements of Mb, we obtain identical constraints on VL m 
and A, but the uncertainties are larger by a factor of ~ 1 .3. 

The band of acceptable values of £!„, and A approximately 
corresponds to a constant growth factor of linear density per- 
turbations between z ~ 0.55 and z ~ 0.05, corresponding to 
the observed factor of 10 evolution of the mass function. For 
small values of fi m , a significant growth factor is achieved 
only if A is large (Carrol et al. 1992). For A = the required 
growth factor is achieved for il m ps 0.4. 

As has become a common practice, we compare our 
(r2 m , A) constraints with those provided by the observations 
of distant type la supernovae and the CMB fluctuations. The 
most robust constraint from the current CMB observations 
follows from the location of the first Doppler peak, which in- 
dicates that the Universe is nearly flat, £! to t = O m + A = 1±0.1 
(with a 0.4 < h < 0.9 prior; Sievers et al. 2002, Rubino- 
Martin et al. 2002). The distant supernovae data require a 
positive cosmological constant. We combine the SN la con- 
straints from Riess et al. (1998) and Perlmutter et al. (1999): 
A- 1.33 Q m = 0.37 ±0. 12. The corresponding 68% confidence 
intervals are shown in the right panel of Fig. [2] Remarkably, 
the constraints from all three methods intersect in a small re- 
gion, n m = 0.27 ± 0.085, A = 0.73 ± 0. 1 1 . 

7. DISCUSSION AND CONCLUSIONS 

We have presented a new method of applying the models 
of the cluster evolution to observations at high redshifts. This 
method is based on easily observable baryon mass function. 
Unlike the studies based on the cluster temperature or lumi- 
nosity functions, our method does not suffer from large un- 
certainties in the M tot -T and M tot -L relations. It uses only 
a physically motivated assumption that the baryon fraction in 
the cluster mass within the large radii is close to the univer- 
sal value, £li,/£l m . Systematic uncertainties in this method are 
smaller and very different from those using the temperature 



or luminosity functions. Until improved observational tech- 
niques are able to normalize the scaling relations to a 10% 
accuracy at all redshifts, modeling of the baryon mass func- 
tion appears to be a promising method for deriving the cluster- 
based cosmological constraints. 

The application of this method to the subsample of distant 
clusters from the 160 deg 2 survey results in a tight constraint 
on the combination of cosmological parameters £!,„ and A. 
This constraint intersects the current CMB and SN la con- 
straints in a small area near il,„ = 0.3 and A = 0.7, showing a 
remarkable agreement of totally independent methods. 

We note an excellent agreement of our cluster-derived cos- 
mological parameter il m with the results based on the X-ray 
measurements of the total mass and baryon fraction in the 
intermediate-redshift clusters (Allen et al. 2002). Also, there 
is a good agreement with those studies of the XTF evolution 
in the EMSS sample which use the normalization of the M—T 
relation based on the X-ray mass measurements (e.g., Henry 
2000). Such consistency indirectly confirms the validity of 
the X-ray cluster mass measurements. 

Our method relies on the assumption that the baryon frac- 
tion in massive clusters does not evolve with redshift. This 
assumption is well-supported by numerical simulations. One 
issue of concern, however, is that these same simulations do 
not reproduce the observed evolution of the Mb - T relation. 
At a fixed T, we observe Mb oc (1 5±0 4 for n m = 0.3 and 
A = 0.7 (Vikhlinin et al. 2002), which is significantly different 
from the expected evolution (M b ,M tot ) oc T 3 / 2 (l +z)" 3 / 2 (e.g., 
Bryan & Norman 1998). This may indicate that the baryon 
fraction in clusters evolves with redshift. However, it is much 
more likely that simulations simply do not reproduce the dis- 
tribution of gas in the cluster central regions. These regions 
usually dominate in the observed emission-weighted T but at 
the same time the gas there is especially prone to processes 
such as radiative cooling which can easily change T but are 
unlikely to modify Mb or M tot . 

The accuracy of our cosmological constraints is presently 
limited by the low number of the distant clusters. When 
new distant cluster samples from the Chandra, XMM, and ex- 
tended ROSAT surveys become available, and more distant 
clusters are observed with adequate exposures by Chandra 
and XMM, even tighter cosmological constraints will be pos- 
sible. Eventually, different methods (cluster evolution, CMB, 
SN la etc.) will start to disagree, which will lead to better un- 
derstanding of the systematics and warrant more complicated 
theoretical models, e.g. those involving non-standard equa- 
tions of state for the dark energy. 

We thank A. Jenkins for providing the software for com- 
puting the model mass functions. Useful discussions with 
M. Markevitch and A. Kravtsov are gratefully acknowledged. 
A. Voevodkin thanks SAO for hospitality during the course 
of this research. This work was supported by NASA grant 
NAG5-9217 and contract NAS8-39073. HQ is partially sup- 
ported by FONDAP Centra de Astrofisica (Chile). 
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FIG. 3. — Left: 68% and 90% confidence intervals on Q„, and A from evolution of the baryon mass function at z > 0.4. Right: Comparison of the cluster 
evolution constraints with those from the distant supernovae and CMB. For each technique, we show only 68% confidence regions. 
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APPENDIX 

MAXIMUM LIKELIHOOD FITTING OF THE OBSERVED MASS FUNCTIONS 

The uncertainty intervals on the low- and high-redshift mass functions were obtained by the Maximum Likelihood analysis of 
the unbinned mass measurements. Our technique is presented in VVPb and briefly summarized below. 



COSMOLOGICAL PARAMETERS FROM CLUSTER EVOLUTION 
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We assume that the observed mass function can be adequately described by the Schechter model, 

f(M) = dN /dM = AM~ a exp(-M /M*) , (Al) 

in the narrow mass ranges of interest. The variation of parameters A, a, and M* within the 3 -parameter N% confidence intervals 
produces a band of the cumulative mass function models which we associate with the N% confidence band for the observed 
cumulative mass function. 

The confidence intervals for parameters A, a, and M* are obtained from the Maximum Likelihood modeling of the unbinned 
observed mass measurements. In the case of negligible mass measurement errors, the likelihood function is (cf. Cash 1979): 

logL= ^ log (/(M,)V(M,)) - f f(M)V(M)dM, (A2) 

where M, are the individual mass measurements, V(M) is the sample volume as a function of mass (computed as described in 
§ 13. H . and the integral is over the entire mass range. The quantity C = -21ogL is statistically equivalent to x 2 (Cash 1979), 
therefore the uncertainties of the estimated parameters can be derived from the standard A\ 2 test. For example, a combined 68% 
confidence region for parameters A, a, and M» is defined by the deviation AC = 3.5 from the minimum (the 68% probability 
point for the x 2 distribution with 3 degrees of freedom). 

In the case of the finite mass measurement errors, the product f(M)V (M) in \A2\ should be convolved with the distribution of 
the measurement scatter. We assume that the scatter follows the log-normal distribution, 

dP(logM) exp(-(logM-log^) 2 /2cr 2 ) 



(A3) 



and the likelihood function is 



dlogM (2w) l / 2 a 
where fj, is the true mass. The convolution of f(M)V(M) with this distribution is 

/•exp(-(logM-logu) 2 /2er 2 ) 
F(M,a)= j (27r)i/V WM^logM (A4) 

logL = Y>g(F (M ; , a { ))- [ F(M,a(M))dM, (A5) 

j JM 

where a-, is the log scatter of the i-th measurement and <j(M) is the trend of the typical measurement scatter with mass. 

The 68% confidence intervals for the cumulative mass function obtained with this technique are shown by grey bands in Fig. [2] 
For the low-redshift mass function, the mass measurement errors are small and we recover the pure Poisson scatter. The derived 
confidence interval for the distant mass function is wider than the Poisson scatter because the masses are less accurate. 

The described technique automatically corrects any biases in the mass function determination which could arise from the mea- 
surement uncertainties (e.g., the Eddington bias). The close agreement between the model (not convolved with the measurement 
errors) and the data in Fig.|2]shows that such biases are small in our case. A bias could be expected because the mass function is 
steep, but is absent because of the rather flat distribution of the observed number of clusters in our sample per log interval of mass 
(this is the product of the intrinsic mass function and the selection function). The effect of the mass measurement uncertainties 
on the mass function is to smooth the distribution of dN /dlogM. If this distribution is flat, the net effect of the measurement 
uncertainties is zero. We verified the absence of biases in our case using Monte-Carlo simulations. 



